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Abstract 

We discuss flavor-violating constraints and consequently possible charged Higgs boson phenomenology emerging from a four- 
zero Yukawa texture embedded within the Type-Ill 2-Higgs Doublet Model (2HDM-111). Firstly, we show in detail how we can 
obtain several kinds of 2HDMs when some parameters in the Yukawa texture are absent. Secondly, we present a comprehensive 
study of the main B-physics constraints on such parameters induced by flavor-changing processes, in particular on the off- 
diagonal terms of such a texture: i.e., from /i — e universality in r decays, several leptonic B-decays {B — >■ ri/, D ^ and 
Ds — >■ Iv), the semi-leptonic transition B — >■ Dtu, plus B X^'y, including B° — B" mixing and the radiative decay Z — ^ bb. 
Thirdly, having selected the surviving 2HDM-111 parameter space, we show that the H~cb coupling can be very large over 
sizable expanses of it, in fact, a very different situation with respect to 2HDMs with a flavor discrete symmetry (i.e., Z2) and 
very similar to the case of the Aligned-2HDM (A2HDM) as well as of models with three or more Higgs doublets. Fourthly, we 
study in detail the ensuing phenomenology at the Large Hadron Collider (LHC), chiefly the cb — ^ production mode 
and the — >■ cb decay channel while assuming t'^Ut decays in the former and t — >■ bH'^ production in the latter, showing 
that significant scope exists in both cases. 
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I. INTRODUCTION 



The main problem in flavor physics Beyond the Standard Model (BSM) [l| is to control the presence of Flavor 
Changing Neutral Currents (FCNCs) that have been observed to be highly suppressed by a variety of experiments. 
Almost all BSM scenarios that describe physics in energy regions higher than the Electro- Weak (EW) scale have 
contributions with FCNCs at tree level, unless some symmetry is introduced in the scalar sector to suppress them. 
One of the most important extensions of the SM is the 2-Higgs Doublet Model (2HDM) 0-3, due to its wide 
variety of dynamical features and the fact that it can represent a low-energy limit of general models like the Minimal 
Supersymmetric Standard Model (MSSM). There are several realizations of the 2HDM, called Type I, II, X and Y 
(acronymed as 2HDM-I 0, 2HDM-II Q, 2HDM-X and 2HDM-Y or inert Types, wherein (part of) the scalar 

particle content does not acquire a Vacuum Expectation Value (VEV) [lTI - [l3 |. In the most general version of a 
2IIDM, the fermionic couplings of the neutral scalars are non-diagonal in flavor and, therefore, generate unwanted 
FCNC phenomena. Different ways to suppress FCNCs have been developed, giving rise to a variety of specific 
implementations of the 2HDM. The simplest and most common approach is to impose a Z2 symmetry forbidding all 
non-diagonal terms in flavor space in the Lagrangian [l5l | . Depending on the charge assignments under this symmetry, 
the model is called Type I, II, X and Y or inert. There are other suggestions for the most general 2IIDM: (i) the 
alignment in flavor space of the Yukawa couplings of the two scalar doublets, which guarantees the absence of tree-level 
FCNC interactions [3 IT7| ; (i i ) th e Lepton Flavor Violating (LFV) terms introduced as a deviation from the Model 
II Yukawa interactions in [I^, [T^; (iii) the 2IIDM-III with a particular Yukawa 'texture', forcing the non-diagonal 
Yukawa couplings to be proportional to the geometric mean of the two fermion masses, gij oc ^mimjXij j20l - l24l )H: 
(iv) recently, a Partially Aligned 2IIDM (PA2HDM) was presented and a four-zero texture is employed therein too, 
so that this newly suggested scenario includes (i) and (iii) as particular cases [27| . 

Therefore, the mechanism through which the FCNCs are controlled defines the actual version of the model and the 
consequently different phenomenology that can be contrasted with experiment. In particular, we focus here on the 
version where the Yukawa couplings depend on the hierarchy of masses. This version is the one where the mass matrix 
has a four-zero texture form [25l . l26j . This matrix is based on the phenomenological observation that the off-diagonal 
elements must be small in order to dim the interactions that violate flavor, as experimental results show. Although 
the phenomenology of Yukawa couplings constrains the hierarchy of the mass matrix entries, it is not enough to 
determine the strength of the interaction with scalars. Another assumption on the Yukawa matrix is related to the 
additional Higgs doublet. In versions I and II a discrete symmetry is introduced on the Higgs doublets, fulfilled by 
the scalar potential, that leads to the vanishing of most of the free parameters. However, version III, having a richer 
phenomenology, requires a slightly more general scheme. Interesting phenomenological implications of 2HDMs with 
a four-zero texture for the charged Higgs boson sector [22|, [2^, [2I] and neutral Higgs boson sector [2^ [s^ have been 
studied. In these works one estimated the order of the parameters Xij — C'(l), including the off-diagonal terms of the 
Yukawa texture. In contrast, a complete and detailed analysis that includes off-diagonal terms of the Yukawa texture 
in presence of the recent data of processes at low energy has been omitted in previous works. Therefore, in this paper, 
we present a comprehensive study of the main fiavor constraints on the parameters that come from a four-zero Yukawa 
texture considering the off-diagonal terms and present their relevance for charged Higgs boson phenomenology at the 
Large Hadron Collider (LHC). The presence of a charged scalar is in fact one of the most distinctive features of 
a two-Higgs doublet extended scalar sector. In the following, we analyze its phenomenological impact in low-energy 
flavor-changing processes within the 2HDM-III with a four-zero texture and constrain the complex parameters Xij 
therein using present data on different leptonic, semi-leptonic/semi-hadronic and hadronic decays. 

If m.H± < nrit — nib, such particles would most copiously (though not exclusively (3ll.l32j) be produced in the decays 
of top quarks via t — )■ H^b [33|. Searches in this channel have been performed by the Tevatron experiments, assuming 
the decay modes — !• cs and tv [13, HH. Since no signal has been observed, constraints arc obtained on the 

parameter space [mjj ± , tan /3] , where tan/3 = v^jvi (i.e., the ratio of the VEVs of the two Higgs doublets). Searches 
in these channels have now been carried out also at the LHC: for — >■ cs with 0.035 fb~^ by ATLAS [3g and for 
TV with 4.8 fb-i by ATLAS [13 plus with 1 fb"! by CMS [1^. These are the first searches for states 
at this collider. The constraints on [m//±,tan/3] from the LHC searches for t — H^h are now more restrictive than 
those obtained from the corresponding Tevatron searches. 

The phenomenology of states in models with three or more Higgs doublets, called Multi-Higgs Doublet Models 
(MHDMs), was first studied comprehensively in with an emphasis on the constraints from low-energy processes 
(e.g., the decays of mesons). Although the phenomenology of bosons at high-energy colliders in MHDMs and 
2HDMs has many similarities, the possibility of TO/f± < mt — together with an enhanced Branching Ratio (BR) 



^ It is well known that, through Yukawa textures [251. 1261 , it is possible to build a matrix that preserves the expected Yukawa couplings 
that depend on the fermion masses. From a phenomenological point of view, the Cheng-Sher ansatz [20II has been very useful to describe 
the phenomenological content of the corresponding Yukawa matrix and the salient features of the hierarchy of quark masses. 
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for — )■ cb would be a distinctive feature of MHDMs. This scenario, which was first mentioned in and studied 
in more detail originally in IH, li^l and most recently in [4l| , is of immediate interest for the ongoing searches 
for t H^b with — >■ cs by the LHC [s^. Although the current limits on — > cs can also be applied to the 
decay — >■ cb (as discussed in in the context of the Tevatron searches), a further improvement in sensitivity to 
t — > H^b with — > cb could be obtained by tagging the b quark which originates from decays [2^, IH, l4lll4^. 
Large values of BR(iy^ — >■ cb) are also possible in certain 2HDMs, such as the "flipped 2HDM" with Natural Flavor 
Conservation (NFC) 0, [13,113 ■ However, in this model one would generally expect m//± ^ mt, due to the constraint 
from b ^ sj (to//± > 295 GeV H^l) so that t -> H^b with — > cb would not proceed. However, in our version 
of the 2HDM-III there are additional new physics contributions which enter b sj, thus weakening the constraint 
on mH± [i^l- We will estimate the increase in sensitivity to BR{H^ — !• cb) and to the fermionic couplings of 
in the 2HDM-III scenario. Further, always in the latter, we will re-visit the possibility of direct production from 
cfo-fusion, where the on-shell state eventually decays to rz^,- pairs, its only resolvable signature in the context of 
fully hadronic machines. 

We now proceed as follows. The formulation of the general 2HDM with a four-zero texture for the Yukawa matrix 
is recalled in section |TT1 The phenomenological consequences of having a charged Higgs field are analyzed in the 
next section in processes at low energy, extracting the corresponding constraints on the aforementioned new physics 
parameters Xij^ by discussing the constraints derived from tree- level leptonic and semi-leptonic/semi-hadronic decays, 
while in section IV we discuss the ensuing light charged Higgs phenomenology at the LHC. Finally, we elaborate our 
conclusions in section V. 



II. THE YUKAWA SECTOR OF THE 2HDM-III WITH A FOUR-ZERO YUKAWA TEXTURE 

In this section, we will discuss the main characteristics of the general Higgs potential and the using of a specific four- 
zero texture in the Yukawa matrices within the 2HDM-III. In this connection, notice that, when a fiavor symmetry 
in the Yukawa sector is implemented, discrete symmetries in the Higgs potential are not needed, so that the most 
general Higgs potential must be introduced. 

A. The general Higgs potential in the 2HDM-III 

The 2HDM includes two Higgs scalar doublets of hypercharge -1-1: $] = (0j", 0°*) and = {4'2 ^(f>2*)- The most 
general SU(2)l x U{1)y invariant scalar potential can be written as [i^ 

l/($i, $2) - M?(*l'Si) + A^2(*2*2) - (m?2(*I*2) + H.c.) + iAi($I$i)2 + iA2($^$2)' + A3($I$l)($J $2) 

+ A4($l$2)($5$l) + Qa5($I$2)' + (A6($I$i) -I- A7($5$2)) ($I$2) + H.C.^ , (1) 

where all parameters are assumed to be reaQ- Regularly, in the 2HDM Type I and II the terms proportional to Ag 
and A7 are absent, because the discrete symmetry $i — > and $2 ~^ —$2 is imposed in order to avoid dangerous 
FCNC effects. However, in our model, where mass matrices with a four-zero texture are considered, as intimated, it 
is not necessary to implement the above discrete symmetry. Thence, one must keep the terms proportional to Ag and 
A7. These parameters play an important role in one- loop processes, where self- interactions of Higgs bosons could be 
relevant [30|. Besides, the parameters Ag and A7 are essential to obtain the decoupling limit of the model in which only 
one CP-even scalar is light, as hinted by current Tevatron and LHC data While these terms exist, there are two 
independent energy scales, v and A2HDM (the scale at which additional BSM physics is required to control persisting 
divergences in the Higgs masses and self-couplings) , and the spectrum of Higgs boson masses is such that m^o is of 
order v whilst 111^0 , m^o and mj/± are all of the order of A2HDM [HI • Then, the heavy Higgs bosons decouple in the 
limit A2HDM ^ according to the decoupling theorem [^. Conversely, when the scalar potential does respect the 
discrete symmetry, it is impossible to have two independent energy scales [Zsj . This implies that all of the physical 
scalar masses lie at the EW scale v. Being that v is already fixed by experiment though, a very heavy Higgs boson 
can only arise by means of a large dimensionless coupling constant A^. In this case, the decoupling theorem is not 
valid, thus opening the possibility for the appearance of non-decoupling effects. Moreover, since the scalar potential 



^ The tJ.i2y -^5' "^6 3.Tid A7 parameters are complex in general, but we will assume that they are real for simplicity. 
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contains some terms that violate the SU{2) custodial symmetry, non-decoupling effects can arise in one-loop induced 
Higgs boson couplings j47l |. 

The scalar potential ([1]) has been diagonalized to generate the mass-eigenstates fields. The charged components of 
the doublets lead to a physical charged Higgs boson and the pseudo-Goldstone boson associated with the W gauge 
field: 

G^. = ^fcp + 4sp, (2) 
= -^tsf)+^fcp, (3) 

with 



„2 2 

m]j± = - -v^iXi + X^+t-'Xe + tpX^), (4) 
spcp 2 1^ 

where we have introduced the short-hand notations, = tan/3, = sin/3 and = cos/3. Besides, the imaginary 
part of the neutral components defines the neutral CP-odd state and the pseudo-Goldstone boson associated with 
the Z gauge boson. The corresponding rotation is given by: 

Gz = cl>licp + 4>lisp. (5) 
^° = -0?/S/3 + 0O,c^, (6) 



where 



TO^o = m^± + iw^(A4 - A5). (7) 

Finally, the real part of the neutral components of the doublets defines the CP-even Higgs bosons and . 
The mass matrix is given by: 

A/«e = i "^'^ ] , (8) 

\^ mi2 m22 y ' ' 

where 

mil = m\sp + w^(Aic| -I- s^As -I- 2s^c^A6), (9) 

m22 = m\cl + v'^{X2s'p + c^As + 2si3CpXr), (10) 

mi2 = -m\spci3 + (^(Aa + Xi)si3Cp + Aec^ + Ays^^ . (11) 

The physical CP-even states, /i° and are written as 

H° = ^%Co, + cf>°2RSa, (12) 

h° = -4>IrSc. + (/-Lca, (13) 



where 



and 



tan 2a ^ ^"^^^ , (14) 
mil - m22 



mnoy: = ^ ( mil + ^22 ± \/(mii - 77122)^ + 4m^2 ) • (1^) 



B. The Yukawa sector in the 2HDM-III with a four-zero texture 

We shall follow Refs. [12, 13, where a specific four-zero texture has been implemented for the Yukawa matrices 
within the 2HDM-III. This allows one to express the couplings of the neutral and charged Higgs bosons in terms of 
the fermion masses, Cabibbo-Kobayashi-Maskawa (CKM) mixing angles and certain dimensionless parameters, which 



4 



are to be bounded by current experimental constraints. Thus, in order to derive the interactions of the charged Higgs 
boson, the Yukawa Lagrangian is written as follows: 



Cy - - \^iQl^iUr + Y^Ql^2Ur + YfQL'^idR + Y^^QL<S>2dR + YILl^^Ir + YlLL^2lRj , (16) 

where $1^2 = {4'i2^ 't'l 2)^ refer to the two Higgs doublets, $1^2 ~ icr2^i 21 Ql denotes the left-handed fcrmion doublet, 
ur and cIr are the right-handed fermion singlets and, finally, Y^2 denote the (3 x 3) Yukawa matrices. Similarly, 
one can see the corresponding left-handed fermion doublet Ll, the right-handed fcrmion singlet Ir and the Yukawa 
matrices F/ 2 fo'" Icptons. 

After spontaneous EW Symmetry Breaking (EWSB), one can derive the fermion mass matrices from eq. (|16p . 
namely 

Mf^^{viY/ + V2Y2f), f = u,dj. (17) 

Wc will assume that both Yukawa matrices y/ and y/ have the four-texture form and are Hermitian [2^ [26j . 
Following this convention, the fermions mass matrices have the same form, which can be written as: 

/ Cf \ 

Mf ={ C} Bf Bf \. (18) 
V B} AfJ 

When _Bg — > one recovers the six-texture form. We also consider the hierarchy | yl^ | 3> | -Bq | , | -Bq | , | Cg | , which is 
supported by the observed fermion masses in the SM. 

The mass matrix is diagonalized through the bi-unitary matrices Vl,r, though each Yukawa matrices is not diago- 
nalized by this transformation. The diagonalization is performed in the following way: 

= VfLMfYfR. (19) 

The fact that Mf is Hermitian, under the considerations given above, directly implies that Vf^ = VfR, and the 
mass eigenstates for the fermions are given by 

u = V^u', d Vjd', I = Vih'. (20) 
Then, cq. (fTT]) in this basis takes the form 

Mf = ^{v,Y/+V2Y/), (21) 

where y/ = Vfi^Y/VjR. In order to compare the kind of new physics coming from our Yukawa texture with some 
more traditional 2HDMs (in particular with the 2HDM-H), in previous works [2^ [23l. [28l - [30| . some of us have adopted 
the following re-definitions: 
2HDM-H-like 

= ^^Md-tan^Y^^, 

V cos p 

/2 

1^" = -^^Mu-cotpYI", 

V sm p 

Yl = Yf{d I). (22) 

These re-definitions are convenient because we can get the Higgs-fermion-fermion coupling in the 2HDM-HI as 
52HDM-III = 52HDM-II + ^9^^'^^ where ff2HDM-ii coupling in the 2HDM-H and IS.gf^'^ is the contribution 

of the four-zero texture. If Ag^^"^ we can recover the 2HDM-II. However, these re-definitions are not unique. In 
fact, there are others possibilities since from eq. (PT|) one can reproduce the 2HDM-I, 2HDM-X or 2HDM-Y as we 
can obtain for any version of 2HDM the following relation: 32HDM-iii ~ -^mDM-any + A'g-^^"^. The other possible 
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re-definitions are: 
2HDM-I-like 





^2 


= — : — 7;Md 
V sin p 


- cot /jy/^ 






= 

vsmp 


- cot I3Y^ 




y' 






2HDM-X-like 


\rd 

^2 


= ''„*^ 

V sm p 


- cot /3Y]^ 




vu 
^2 


V sm p 


-cot/3Y]" 




n 


= Yi'^Cd /) 




2HDM-Y-like 


yd 


^ ^ M. 
V cos p 


-tan/gfj'* 




vu 
^2 


V sm p 


- cot PY^ 




^2 







(23) 



(24) 



(25) 

After spontaneous EWSB and including the diagonalizing matrices for quarks and Higgs bosont|f|, the interactions of 
the charged Higgs bosons and neutral Higgs bosons {(fP = h^, , ) with quark pairs for any parametrization 
2HDM-(I,H,X,Y)-like have the following form: 



£fifj<t> 



2^/2Ml 



w 



1=1 

ym„, Su-f{Y) 
Zmi^ 5ij - f(Z) 



(1 + 7^) 



9 

f V2Mw 
\ 9 



(VcKM)i,(l-7') \ d,H 



(26) 



n{Z) 



{l + 'y'')ljH+ + h.c. 



2M« 



rrid-di 



f{X) g \md, ' 



+ 



-x5,, + /(x)^f!!i^) (?;V)).- 



^" ' f{Y) T ) ^ "^^'^'^ 



YS,, + f(Y) 



f{Y) g {m l^"^" 
9 K'Triui 



(C( 



7^" U, 



+mi.li 



+i 



(d + ^Ci/) %/2 / 

(Cl/ - ^CL) 72 fmw] 
f{Z) g \mdj 

-Z5., + nZ)^ f^) (F^,,,)., 
g \mdj 



{Yn{Z))i 



The details of both diagonalizations are presented in Ref. [22 
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where Vqkm denotes the mixing matrices of the quark sector, the functions f(^x) and n{x) are given by: 



2 



fix) = Vl + x- 

, , (2 if a; = tan /3, 
"(") = il ifx = cot?. (27) 



the parameters X,Y,Z are given in Refs. jl. Isl-fiol HgI. l4ll. lisj and the factors are presented in Ref. Following 
this notation we can list the parameters for the framework 2HDM-(I,II,X,Y)-like through Tab. H] 

Following the analysis in one can derive a better approximation for the product Vq Y,'^ , expressing the rotated 



matrix Y^, in the form 



— ixlh = ^^-—ix%^'^' (28) 



where the x's are unknown dimensionless parameters of the model, they come from the election of a specific texture 
of the Yukawa matrices. It is important to mention that eq. (j28p is a consequence of the diagonalization process of 
Yuwaka matrices, assuming the hierarchy among the fcrmion masses (see Ref. [131), namely, the Cheng-Sher ansatz 
is a particular case of this parametrization. Besides, in order to have an acceptable model, the parameters x's could 
be 0(1) but not more, generally. Recently we have calculated the fit of Yukawa matrices including the CKM 
matrix, and wc find that the parameters off-diagonal arc 0(1) (e.g., X23 — 10): therefore we cannot ignore all of 
these [43|. Besides, in Ref. [5Q|, they study the general 2HDMs considering rcnormalization group evolution of the 
Yukawa couplings and the cases when the 2^2-symmetry is broken, called non-diagonal models (e.g., the models with 
a structure incorporating the Cheng-Sher ansatz). It is interesting to note that it is actually the ofi^-diagonal elements 
in the down-sector that become large whereas the ones in the up-sector x''(m) ^ 0.1, assuming the conservative 
criterion ^ 0.1, where fx is the rcnormalization scale. On the other hand, the FCNC processes at low energy are 
going to determine bounds for these parameters with high precision, aspect which is studied in this work. In order 
to perform our phcnomcnological study, we find it convenient to rewrite the Lagrangian given in eq. (|26l) in terms of 
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the coefRcients as follows: 
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2V2Mw 



1=1 



+ TO„, Sil - 

9 



m 

f{Z) 



( fi-^) \ 

(Vckm);; (x ma, Sij - '—^ y/md,md, xfj) (1 + 7^) 
y/mi^nid, x\j ) (1 + l^)ljH+ + h.c. 



2M, 



w 



V2 



^diKh"ij ^ 77^\ 7^ X 



f{X) V2 



f{X) 72 



if" 



-rnd^X5,,+f{X)^l^^X% 



^uiC,h"i] 77^7^; \ ) Xij 



f(.Y) 



^/2 



+i 



-m.,yj,,- + /(y) ^ ^ ' xl 



mii^H^ij 777r\ 7^ — Xt 



f{Z) V2 
rniiihOij H 7777^ 7^ — Xi 



(29) 



(30) 



where we have redefined [xil^^- = Xtj, l-^-zuj -^tj 
UidjH^ and UidjH^ arc given by: 



xfj and [xa],^ = x\j- Then, from eq. (EH), the couphngs 



9H+Uidj 



iff 



2Mh 



2^Mv 



2^Mv 



(31) 
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2HDM-III 


X 


Y 


Z 


?/l 


?h 


pd 


All 
?/f 


cd 
t,H 




2HDM-I-like 


-cot 13 


cotp 


-cotp 


Ca/Si3 






Sa/Sp 


Sa/sp 


Sa/sp 


2HDM-II-like 


tan P 


cotp 


tanp 


Cajsp 


— Sa/C0 


— Sa/Cp 


Sa/Sp 


Ca/Cp 


Ca/cp 


2HDM-X-like 


-cot/3 


cot P 


tan P 


Ca/Sp 


Ca/Sp 


-Sa/Cp 


Sa/sp 


Sa/sp 


Ca/cp 


2HDM-Y-like 


tan P 


cot/3 


-cotp 




-Sajcp 


Cq / Sp 


Sa/sp 


Ca/cp 


Sa/sp 



TABLE I. Parameters X, Y and Z defined in the Yukawa interactions of eq. p6|) for four versions of the 2HDM-III with a 
four-zero texture, which come from eqs. ((22}-((25]). Here Sa = siuQ, Ca = coso?, sp — sin/3 and cp — cos/3. 



where h{j, H{j, A{j, Sij and Pij are defined as: 



Ki = H 7^T777^\lz:—Xip Ki = h,Ad ^ I, X ^ Z), 



V2f{X) y rud^ 



rrd _,d. i^t + X^i) fmj- d 



-XSi, 



V2f{X) y rud. 

f{X) [7^ Al^^AfJd^l,X^Z), 



V2 y md 



V2f{Y) V^«. 
V2fiY) 



'Xiji 



4" — —V/i 

Sij = Xij + rriu^ , Pij = nid^ X^j - m^^ Yij, 



f{Y) /Jf^ 

Xi 



with 



Xij — ^(Vckm)i 
1=1 

y. = t 



'rud, f{X) fW^ 



1=1 



YSa - 



a/2 V 



■X^l 



(VcKM);j- 



For the case of leptons S\j = Pj^ we have 



S\j — mi. Z\^ , 
mii 



'Z'lj - 



mi, 



Sh - 



V2 V 



Then, the couplings i^ijH^ and l^vi.H are given by 



(32) 



(33) 



(34) 



(35) 



(36) 



In order to compare these couphngs with previous works [1, Isl-fiol. 111,111 

, we find it convenient to define the couphngs 
UidjH~^ and UidjH~ in terms of the matrix elements Xij^ Yij and Zij. Following the definitions (l32])-([35]) we obtain 
the following compact expression for the interactions of Higgs bosons with the fermions: 



V2. 



^/2\ 



Cf^f,^ = _ ) Xlji^ {md,X,jPR + m^Y.jPL) d, H+ + 1^!^Z,,WIrH+ + H.c. 



(37) 



When the parameters Xij = 0, '^^ obtain Xn = X22 = -^33 = X (similarly for Y and Z) and one recovers the 
Yukawa interactions given in Refs. j^. [sl-fiol l4l|. Besides, in order to hold consistencies with the MHDM/A2HDM 
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3, we suggest that this Lagrangian could represent a MHDM/A2HDM with additional flavor physics in the Yukawa 
matrices as well as the possibility of FCNCs at tree level. Returning to the Lagrangian given in eq. ([37|) . when the 
parameters xfj are present, one can see that Xu ^ X22 7^ ^ X and the criteria of flavor constraints on X cannot 
be applied directly to Xij (the same for Y and Z), but the analyses for low energy processes are similar. Below we 
shall discuss more about various aspects of the model. Finally, it should be pointed out that parameters X , Y ^ Z , 
and Xij a-re arbitrary complex numbers, opening the possibility of having new sources of CP violation with tree-level 
FCNCs. 

Previously, in Ref. [13] j the flavor constraints of the 2HDM-III with a six-zero texture were studied, finding 
interesting results that we can use. However, we should compare their results and ours so as to distinguish the two 
parametrizations. Firstly, the six-zero texture assumed in [2J| has been disfavored by current data on the CKM 
mixing angles [2^[5l|. Hence, we focus here onto the four-zero texture, which is still acceptable phenomenologically 
and of which we consider the non-diagonal terms of the Yukawa matrices. Secondly, in order to unify notations we 
relate the parameters Xf^ of [IJl with our parameters Xij, Yij and Zij as given in eqs. (p4| and ((35)) . as followfQ: 



ftj = Azj , 



(38) 



(39) 



where pfj was introduced following the Cheng-Sher ansatz, considering AfJ ~ 0(1). If we compare this with cq. p7p . 
after using eqs. (p4)) and ([35]). we obtain the following relations: 



- 



Z . I d,, — Xi, 



mi, 



V2 



(40) 



and 



(41) 



In essence, in the remainder of our work, we assume that our model could represent an effective flavor theory, wherein 
the Higgs fields necessarily participates in the flavor structure and has the same features as those of renormalizable 
flavor models (52| - [55j . In this type of scenarios, a horizontal flavor symmetry, continuous or discrete, is added to the 
SM gauge group symmetry in such a way as to reproduce the observed mass and mixing angle patterns by only using 
renormalizable terms in the Lagrangians. This requirement has two immediate and interesting consequences: firstly, 
there must be more than one SU (2) doublet scalar; secondly, at least some of them must transform non-trivially 
under the flavor symmetry (sgI . [57| . 



III. FLAVOR CONSTRAINTS ON THE 2HDM-III WITH A FOUR-ZERO YUKAWA TEXTURE 

In this section we will analyze the most important FCNC processes that are sensitive to, in particular, charged Higgs 
boson exchange, the primary interest of this paper, as well as effect of a (neutral) SM-like Higgs boson (we assume 
that its mass is m/jO ~ 125 GcV). Starting from measurements obtained from from these processes we constrain the 
new physics parameters xfj that come from four-zero Yukawa texture. Finally, we study the possibility of obtaining 
a light charged Higgs boson compatible with all such measurements. We will address the various experimental limits 
in different subsections. 



We adopt the description of the Yukawa sector presented in Ref. [24 
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FIG. 1. Considering the constraint from /i — e universality in r decays, we show the allowed region (orange color) for X22 a^nd 
X33 when Z takes values of 10 (left-panel), 40 (center-panel) and 80 (right-panel). Here, 90 GeV < mfj± < 130 GeV. 



A. — e universality in r decays 

The r decays into fiPp^Vr and eDcMT produce important constraints onto charged Higgs boson states coupling to 
leptons [ssf . through the requirement of /i — e universality. The consequent limits can be quantified through the 
following relation i63| : 



BR{t — > eDeVr) girni/ra^) 



= 1.0036 ± 0.0020 



where g{x) = 1 - 8.t^ + 8a;^ 
following relation 



12x^ logx. Following [8|, in our case, the request of /i — e universality imposes the 



BR{t ^ ^iD^Ur) f {ml/ml) ^ ^ ^ ^ _ q 25R 
BR{t — > eVgi/r) .f{fnfj,/^r 4 ' 



(42) 



where R is the scalar contribution parametrized through the effective coupUng, see eqs. ((35|) and ((36 



R= \ Z- 



■2^22 — 



33 ^22 



z f^^K^ ' 



(43) 



One can see that R is symmetric in the two parameters X22 and 
obtain the following explicit constraint: 



X33. Following the analysis of Ref. [61[, we can 



\ZiiZ; 



< 0.16 GcV"^ 



(95% CL). 



(44) 



We show in Fig. [T]the constraints on X22 and X33 with Z = 10, 40, 80. One can see that, for small Z values, the 
allowed region for X22 and X33 is large whereas, when Z instead grows, the allowed region for theses parameters is 
smaller in comparison. The constraints becomes most restrictive when Z is large and we have a very light charged 
Higgs boson, between 90 and 130 GeV. The plot also shows that X22 and X33 could be simultaneously 1 and — 1, 
respectively, and the more favorable region is the one where X22 = X33 = 1-5 for 0.5 < Z < 100. When Z is large, if 
X22 = 1; one can see that 0.5 < X33 ^ 2.5 (the same happens when X22 and X33 are interchanged). Further, in Fig. 
[2] we present the plane [mH± , X] and the allowed region is shown for the cases X33 = and X22 = (l^ft panel) and 
X22 = 0.1 and — 20 < X33 ^ 20 (right panel). In the left panel we present the red region (without contributions from 
the parameters of flavor physics |xy |)j which is allowed by /i — e universality in r decays: e.g., for to^± < 120 GeV we 
must have the constraint X < 50. In the right panel we show two regions: here, the blue(blue&gray) one is allowed 
for the cases 5 < IXu KIXijl ^ 5). The blue region is clearly the more restrictive one of the two and could become even 
smaller while the Ixyl's grow. In the case shown, we can get that 77i^± < 150 GeV for X < 20. Conversely, with 
both regions combined, blue&gray, which represent the portion of parameter space allowed for 0.8 < \xij\ ^ 2, we 
see that the model is more favored, because it opens up larger regions in the plane [mH±,X]. One can see this, e.g., 
for X < 80, as the bound for the charged Higgs boson mass is now given by to^± > 100 GeV, that is, not dissimilar 
from the previous case (when X < 20). 
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FIG. 2. Considering the constraint from /i — e universality in r decays, we show the allowed region for the plane mjf± — X: 
left panel represent the case X33 ~ and X22 ~ 0, right panel represents the case X22 = 0.1 and —20 < X33 5- 20 (blue region 
is for 5 < ixiji) blue&gray region is for jx'jl < 5). The same occur when X33 = 0.1 and —20 < X22 1^ 20. 



B. Leptonic meson decays 

The leptonic decay of a charged meson, M li^iAs sensitive to iJ+ exchange due to the helicity suppression of the 
SM amplitude. The total decay width is given by [6ll.[6^: 

r(M,, ^ Ij,) ^ G%miflj\V,,f'^{l + (Se™)|l - A,,|, (45) 

an 

where i, j represent the valence quarks of the meson. Vij is the relevant CKM matrix element, /m is the decay 
constant of the meson M (the normalization of the meson decay constant correspond to /,r = 131 MeV), Sem denotes 
the electromagnetic radiative contributions and A^- is the correction that comes from new physics information. In 
particular, for the 2HDM-III employing a four-zero Yukawa texture, the leptonic decays receive a contribution from 
charged Higgs bosons in the following form; 



A.= f^U..ffe^^), k = 2,3. (46) 



2 

V mH± ) V ("^M, + rndj ) 

In the more general 2HDM-III the Aij correction can be a complex number. As is pointed out in [6l] |. in some Two 
Higgs Doublet Models (2IIDM's) with natural flavor conservation the correction A^- (with x's = 0) is predicted to be 
positive (in 2HDM-I) or negative (in 2HDM-X), while can have either sign in 2IIDM-II and 2IIDM-Y, depending on 
the decaying meson, whereas it is absent in the inert Higgs scenario. 

We focus on decays of heavy pseudoscalar mesons D — > fii/, B tv and Ds — )• ^v, tv, which have been measured. 
In B and D decays the function A^, one can neglect the contribution proportional to the light quark mass because 
Tnu/frib < md/rric 0(10""^). Hence the functions Aij for D — )■ /iiy and B — > tv, respectively, are given by: 



A„ 



ml 


ry 

Vcd 


ml 1 


' Z - 










ml 


ry -'^13 

Vub 


ml 






ml± 


ml± 


^/2 



f{Z) f{Y) l^VtafiY) „\ 

^ut ~ 3^ ^33— - ^ - ^X33 ) ( ^X33) - (47) 

Apparently, the factor ^/rnt/rn^ in Ad could be considered as a dangerous term, which could make the theoretical 
predictions deviate from the experimental results, however, the term Vtd/Vcd reduces this possible effect. Similarly, 
this happens when one wants to fit the four-zero texture of the Yukawa matrices with the CKM matrix. Since 
experimental results of B{D^ iiv), which were measured by CLEO collaboration [g^I, the authors of Ref. [6l| 
found the following constraints at 95% C.L. for any model: 0.8 < |1 — A„f,| < 2 and 0.87 < |1 — Acd\ < 1.12. 
Considering those constraints, we can get the allowed circular bands in the Z 22Y21 / (m^^Vcd) and Z^^Xi^ / (ml^Vub) 
complex planes. Our numerical analysis obtained from the decays B ^ tu and D — > pLV is shown in Fig. [3l which is 
consistent with the results of [Glf when the parameters x's are absent. For instance, we also find the real solutions 
are ^33^i3/(m^±Kb) € [-0.036,0.008] GeV^^ or [0.064,0.108] GeV-^ from the B tv, and Z22Y2i/{ml±Vcd) e 
[-0.037,0.035] GeV-2 or [0.535,0.609] GeV~2 from the D fiv. In Fig. Hwe show the allowed region for the plane 
[X22j X23]' assuming that X22 € [O-l, 1-5], and considering the bounds for D ^v. One can see that X23 ^ [0.75, 1.25] 
when X22 = 1 30 < \Z\ = \Y\. For the cases Z » Y ai Y » Z the permitted region is larger than for 
1 << \Z\ = |y| and 1 < X23- Therefore, 1 < X23 are allowed parameters for the leptonic decay of D mesons and the 
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consequences on the phenomenology of charged Higgs bosons could be an important probe of the flavor structure of 
the Yukawa sector. Otherwise, for the low energy process B riy we can get bounds for the parameters of the Yukawa 
texture pertaining to the d-quark family. In Fig. [SJ we show the allowed regions in the plane [X221 X23] the following 
cases: X >> Z (left panel), Z >> X (center panel) and Z,X >> 1 (right panel), with 80 GeV < m^± < 160 GeV 
and considering 0.1 < X22 < 1-5. We can see that the non-diagonal parameter xis is more constrained than X23- For 
the case Z, X >> 1 and X22 — -^i ^ — — ^0 we found X23 ^ [ — 0.35, — 0.2] or X23 ^ [OjO.2], so that this case 
could correspond to a 2HDM-II-likc scenario (see Tab. H]) when tan (3 is large. This scenario is more constrained when 
X = Z > AO and the bound for IX23I < 0.2 is obtained. Another interesting case is when Z >> X and X22 = 1, for 
X = 0.1 and Z = 80 we obtain the following allowed regions: X23 ^ [^1-8, —1-2] or X23 S [—0.2, 0.6], in this scenario 
it is therefore possible to obtain the constraint IX23I = 1- When X >> Z wc get a wider permitted region for X23J 
defined in the interval (—7, 2). 

From Ds — > iJiv,Tv decays the constraint 0.97 < |1 — ^cs\ < 1-16 is obtained in Ref. [6l| and one can get the real 
solutions t^cslm^B ^ [—0.044,0.008] GeV ~^ or [0.545, 0.598] GeV Here the situation is somewhat more complex, 
because ms/rric ~ 10% and we cannot neglect s-quark effects, so the expression for Acs is: 

A f'^D,\^ry f Y22'mc + X22ms\ „ 

Acs = ^ Zkk -f^T ^ fc = 2, 3 , 48 




mh fix) 

^cb 7=- 



(49) 



With this information, we can establish a correlation among the parameters that come from D fiiy and B — >■ tu. 
Considering the information from B — ^ riy, D — > fiiy and Ds tv, fiv, we show in Fig. [6] the constraints for the 
non-diagonal terms of the Yukawa textures X23 ^-^id X231 assuming 0.1 < X22 = X33 ^ 1-5: as well as X22 — X22 — 1- 
Wc present in Tab. |ll]a set bounds for these parameters in several scenarios, which are shown in Tab. ID Combining 
results from the table, one can derive general constraints for IX23I < 0.15 and IX23I ^ 1-5 for almost all scenarios. 
Only in the 2HDM-Y-like version one can obtain a less stringent bound for X23- 



C. Semileptonic decays B — !> Drv 

Purely leptonic decays of mesons interwine EW and QCD interactions. However, the role of strong interaction 
materializes only in the presence of a decay constant, to be assessed through theoretical methods. Semi-leptonic 
decays are complicated to describe since they involve form factors with a non-trivial dependence on the momentum 
transfer. If the form factors are known with sufficient accuracy, scmi-leptonic BRs start becoming stringent constraints 
on new physics models. The BaBar and Belle experiments published the first measurements of B{B — > Dtv) (63.[65j. 
Recently, using the full data set collected by BaBar, the update of BR(i? Dtv) and BR(i? D*tv) was presented 
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FIG. 4. The most constrained region for X22 s^nd X23 from D ^ fiu for the following cases: Z, Y >> 1 (left), Z >> Y (center) 
and Y » Z (right), with 80 GeV < m„± < 160 GeV. We assume that 0.1 < X22 < 1-5. 
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FIG. 5. The most constrained region for X22 a-nd X23 from B ^ tu for the following cases: X » Z (left), Z » X (center) 
and Z, X»l (right), with 80 GeV < ?n^± < 160 GeV. We assume that 0.1 < X22 < 1-5. 



in [6g, from where it is clear that the 2HDM-II is disfavored. Since this model cannot explain R{D) and R{D*) 
simultaneously (and for B ^ tv & high fine tuning is needed), where R{D*) are the ratios 

R{D*) = BR(.B ^ D*Tv) /BR(S ^ D*lv) (50) 

with 

i?(i:>) = 0.44 ± 0.058 ± 0.042, (51) 
R{D*) = 0.332 ± 0.024 ± 0.018. 

However, lately, in Ref. [g^l, it was shown that one can simultaneously explain R{D) and R{D*) in the 2HDM-III 
with a general flavor structure, where the non-diagonal terms from the u-quark sector are relevant. 

Following the analysis of Ref. [i^l, an interesting observable is the normalized BR, Rb^Dtu = BR(i? — 
DTi')/BR(i? — > Dev), which corresponds to a & — )- c transition, with a CKM factor much larger than the purely 
leptonic B decay. One can write this term as a second order polynomial in the charged Higgs boson coupling to 
fcrmions, as 

Rb^Dtv = ao + al("^B - m^jhs + a2(™| - m^ofSl^, (52) 
where the factor ^23 is determined by the coupling H^Uidi, where the general expression for 5ij is given by 



- -"■'■■■-^ ). (53) 



The polynomial coefficients in eq. ((52|) are given in Ref. |62| as: 



ao = 0.2970 + 0.1286dp^ + 0.7379dA, 

ai = 0.1065 + 0.0546dp2 + 0.4631dA, (54) 
a2 = 0.0178 + O.OOlOdp^ + o.0077dA, 
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X'i^ Xn Xn ^^23 

FIG. 6. The most constrained region for X23 vs. X23 from B — >■ ru, D ^ ^iv and Da — >■ Iv for the following cases: \Z\ — \X\ = 
|y| >> 1 (left), Z » X,Y (center-left), X »Y,Z (center-right), and \Z\ = \X\ » Y (right), with 80 GeV < m^± < 160 
GeV. We assume that 0.1 < X22 < 1-5 and X22^ = — 1- 




FIG. 7. The most constrained region for xis vs. X23 from B ^ tv and B Dtv for the following cases: \Z\ = |X| = |y| >> 1 
(left), Z » X,Y (center-left), X »Y,Z (center-right), and \Z\ = |X| >> Y (right), with 80 GeV < m„± < 160 GeV. We 
assume that 0.1 < X22 1-5 and x^i'' ~ X33'* — 1- 



where dp^ — — 1.18 and dA = A — 0.046 are the variations of the semi-leptonic form factors p^ and A [gI, [69l |. 
Similarly to the leptonic process Ds Iv, we can establish a correlation among the parameters that come from 
B — >■ Dli' and B — >■ riy. In all cases, we consider simultaneously R{D) and R{D*). One can then constraint the 
non-diagonal terms of the Yukawa texture, X23 ^^'^ X23' by assuming 0.1 < X22 = X33 — 1-5 as well as X22 = X22 = 1- 
We can then show in Tab. |n]a set of bounds for these parameters in several scenarios. By combining results from this 
table, one can derive general constraints for IX23I ^ 0-15 and IX23I ^ 1-5 for several scenarios (see also Fig. [7]). Again, 
the 2HDM-Y-like version cannot offer a bound for X23 easily. One can see that the 2HDM-III with a Yukawa texture 
can avoid the constraints of the factor Rb^Dtd and can thus appear as rather exotic physics, i.e., very different from 
the traditional 2HDMs with NFC. In particular, decay channels involving cb could be relevant and be searched 

for in the transition t — > H^b [H, |4l|, if the state is sufBciently light. 



2HDM-III versions 


xUB^tu) 


xhiDs^liy) 


xhiB-^Dru) 


X23 (combination) 


{X,Y,Z) 


2HDM-I-hke 
2HDM-II-like 
2HDM-X-like 
2HDM-Y-like 


(-0.35,-0.15) or (0,0.15) 
(-0.35,-0.2) or (0,0.2) 
(-7,2) 
(-1.8,-1.2) or (-0.2,0.6) 


(-1.5,0.9) 

(-2,27) 
(-4, 14) 
(-40,50) 


(-0.05,0.45) 
(-9.6,1.2) 

(-3.8,0.47) 
(-16,50) 


(-0.05,0.45) 
(-2,1.2) 

(-4.2,0.47) 
(-16,50) 


(20, 20, 20) 
(20, 0.1, 20) 
(0.1, 0.1, 50) 
(50, 0.5, .5) 



TABLE II. Constraints from B — Drf, Ds — > tu, pu and B ^ tu decays. We show the allowed intervals for X23'' constrained 
by each low energy process, according to the different scenarios presented in Tab. Uas well as the combination of constraints 
for the X23 parameter. We assume 0.1 < X22 ~ X33 1-5 as well as X22 = X22 = 1- Taking 80 GeV < rnjj± < 160 GeV and 
specific values for the X, Y and Z parameters given in Tab. HI 



D. B ^ As 7 decays 

The radiative decay B — )• has been calculated at Next-to-Next-to-Leading Order (NNLO) in the SM, leading 
to the prediction BR{B Xs-/)sm = (3.15 ± 0.23) x 10"'' Hil. In the 2HDM the decay amplitude is known at NLO 
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[illlSlzii and in the 2HDM with Minimal Flavor Violation (MFV) [t^ while, only very recently, NNLO results have 
been presented for both a Ty pc- I and Typc-II 2HDM [tII . The current average of the measurements by CLEO , 
Belle [iIItI, and BaBar (HH] reads BR{B X^)\e^>i.6 goV = (3.37 ± 0.23) x IQ-'*. 

In this subsection we show the constraints on the ofF-diagonal terms of the four-zero Yukawa texture of the 2HDM- 
III through a general study of the processes B — >■ X^^. We first start with a digression on Wilson coefficients entering 
the higher order calculations. 



1. NLO Wilson coefficients at the scale 

To the first order in a^, the effective Wilson coefficients at the scale = 0{Miy) can be written as (isl [80| 

C^M = C^'^^'ipw) + ^^G'r^M . (55) 



' 47r ' 

The LO contribution of 2HDM-III to the relevant Wilson coefficients at the matching energy scale take the form 

MM, 



^33^32 



VtbVts 



^i7,S).YYiyt) + 



^33^32 



VtbVts 



(56) 



where yt = m?/TO^±, SC°^^^^ (nw) = C°^%HlJ-w) - 8),SAf(MvK), and the coefficients C^^^s^snil^w)^ C'[l,_g) (yt), 



C"y xY^yt) ^^^^ known, which are given in [43|, |80| and 



^33^3*2 



VtbVts 

^33^3*2 
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mt \ VtbJ 



X23 
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m 

V2 



X33 



VtsJ V2 



X23 



VtbVts 



X 



^li^^A^ [^(Vts\f{X) . 



V2 



■X33 



mb \ Vtb J V2 



■X23 



^X33 



mt \ VtsJ V2 



X23 



The NLO Wilson coefficients at the matching scale in the 2IIDM-III can be written as [43 

15 + 6 In- 

III. 

Eh 



(7,8) 'yf^W) 







^33^32 



VtbVts 



lAv 








1 + 


^^33^32 




VtbVts 


(* = 2, 3,5,6), 


^(V,8),yF(/^w') + 



^33^32 



VtbVts 



^(7..S),XY{^^w)^ 



(57) 
(58) 

(59) 

(60) 

(61) 
(62) 



where the functions on the right-hand side of eqs. ([50)1 and are given in Ref. [4^ jSO]. The contributions of our 
version 2IIDM-III to the B — decay are described by the functions C°^^(/^iv) (* = 7, 8 and j = (YY, XY)) and 
the couphngs Fg'^, ¥33* and X33. 

The contribution from new physics to i? — > Xs^ also depend on the magnitude and sign of the couplings X and 



Y (see Tab. For instance, for the conventional Type-II 2IIDM, we have \Y\ 



2IIDM-III under study here we have 

5^331^3*2 



, XY* 



VtbVts 

-^^33^3*2 



VtbVts 



\t 



A 



rncfVcb\f(Y) 



X23 



mt \ VtbJ 
l^(Vts\f{X) / 
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mb \ Vtb J \f2 



v,s\![y) 

mt \ Vts) \/2 

^c(v,s\!{y) 



X23 



mt \ Vts) \/2 



-X23 



1, whereas for the 

(63) 
(64) 



where and A^, expressed in eg. (PH)) . are parameters defined in a version of the 2IIDM-III without off-diagonal 
terms in the Yukawa texture [l^, Isolpl. Again, when the off-diagonal terms of the four-zero texture of Yukawa matrices 
are absent, we recover the results mentioned. When the heavy charged Higgs bosons is integrated out at the scale /^vi^. 



5 In the version 2HDM-III of [ii,!!^,!!!, one has = X^b, X^t = ^tt- 
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the QCD running of the the Wilson coefficients Ci{fiw) down to the lower energy scale /ifc = 0{mb), after including 
the new physics contributions, the same happens in the SM. Thence, for a complete NLO analysis of the radiative 
decay B — > Xgj only the Wilson coefficient C^^^ {iib) has to be known, which is: 

C;/^(m.) - C^f'ip,) + , (65) 

where the functions C^' '^■^ {iJb) and C]''^^^ {^ih) as functions of C'^ ^{^w) and their complete expressions are given in 



E. BR(B Xs-i) 



The BR of the inclusive radiative decay B — > Xg^ is given by: 



BR(B ^ Xsi)lo = BsL 



Vcb 



TTe{z) 



at the LO level and 



NLO 



^B, 



SL 



Vc 



cb 



6a„ 



t:9{z)k{z) 



[\D\^ + A + A\ 



(66) 



(67) 



at the NLO level, where Bsl = (10.74zb0.16)% is the measured semi-leptonic BR of the B meson [63|, aem = 1/137.036 



is the fine-structure constant, z 



pole 



is the ratio of the quark pole masses, 0{z) and k{z) denote the phase 



space factor and the QCD correction [8lJ for the semi-leptonic B decay and are given in [43|, |80|- The term D in eq. 
(|67l) corresponds to the sub-processes 6 — > 57 [ij 



(68) 



where the NLO Wilson coefficient C^^{^b) has been given in eq. (|65l). and the function V(ub) is given by [H, |80|- In 
eq. (j67p . term A is the the correction coming from the bremsstrahlung process b — sjg j82l |. Now we are ready to 
present numerical results of the BRs in the 2HDM-III. We employ the central value of the input parameters given in 
Refs. [i^ls^l- For the values of the matching scale and low energy scale, we take fiw = My/ and mb/2 < /it < 2mf,. 



Following the recent analysis of Refs. [6l|, |83[ and using standard values [43|, |80j for the charged Higgs boson mass 



(80 GeV < mH± < 300 GeV), we can establish the following constraints: 



^33^3*2 



VtbVt 



ts 



< 0.25, 



-1.7 < Re 



-^^33^3*2 



VtbVt 



ts 



< 0.7. 



(69) 



Since 



VtbVts 



< 0.25, we show in Fig. [S]the allowed area in the plane X33 — X23' ^'^^ ^^'^ cases Y << 1 (left panel), 
y = 1 (center panel) and y = 10 (right panel). One can then extract the bounds —0.75 < X23 — ^0.15 for X33 = 1 and 



0.4 < X23 ^ 0.9 for X33 = —1; both when Y << 1. Otherwise, using the second constraint —1.7 < Re 



VtbVts 



< 0.7, 



we can obtain the interval permitted for X237 assuming the allowed interval for X23 from B tv and X33 = 1 = 
X33 = 1. In Fig. [S]one can get the allowed area for some scenarios of Tab. |IT1 We can, e.g., obtain X23 G (—0.55, —0.48) 
for the case X = 2Q and y = 0.1 (left panel). An interesting scenario for the 2IIDM-III is the 2HDM-X-like one, 
where the allowed region is larger than in other scenarios, with X23 ^ (—2.2,0.45) and X23 ^ (~'''j~2) (using the 
constraint coming from B — > tv), so that one can avoid the most restrictive constraints hitherto considered. 



F. B° - B° mixing 

In generic 2HDMs, the charged Higgs boson contributes to the mass splitting AMb^- Previously, in Ref. [sj], the 
new physics contribution to AMb^ was calculated in the 2HDM-III at the LO and was presented the constraints on 
the Xtt — rnB± plane, by using the measured Xd = AMb^/Tb value. However, in that work did not include the effects 

of the large uncertainty affecting the non-perturbative parameter Jb^ Bb^ and the new physics contribution to the 
parameter rjB [soj . as well as the possibility of non-diagonal terms in the Yukawa texture. Remembering that in 
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FIG. 8. The allowed region for xxi ^s. X23 from B — > (using the constraint 



Y33Y, 



22. 



VtbVts 



< 0.25) and B° - mixing 



(considering the constraint 
80 GeV < m„± < 200 GeV 



VtbVtd 



< 0.25) for the following cases: \Y\ « 1 (left), Y = 1 (center) and F = 10 (right), with 



-0.7 -0.6 -0.5 -0.4 -0.3 

Xii 

A 



0.22 0.23 0.24 0.25 0.26 0.27 0.28 
Xl3 




FIG. 9. The allowed region for X23 vs. X23 from B Xs^ (using the constraint — 1.7 < Re 



1 

VtiVts 



< 0.7) for the following 



cases: X = 20 and F = 0.1 (left), X = F = 20 (center) and X = F = 0.1 (right), with 80 GeV < m^± < 200 GeV. We assume 
X33 = 1, X33 = 1. 



Ref. [63 the general flavor structure for the 2HDM-III has been found more consistent with various other B physics 
constraints, we are motivated to test this model assumption also against limits coming from — B^ mixing. 

Being that AM^^ has been measured with very high precision (68j . we utilize this quantity directly instead of 
the parameter Xd- We will compute for our version of 2HDM-III, the charged Higgs boson contribution to the mass 
difference AMb^ at the NLO level, considering the calculations of Refs. [lo, HE]- In the 2HDMs, the NLO the mass 
splitting AMb^ is given by (HHi] 



AMb, = ^mBM^\VtdnVtb\'iBBj'BjVB{xt,yt)S2HDMi 



xt,yt) 



where Xt = m'j{Mw)/M^^, yt = 'm1[Mw)/m?u± and 



as{Mw) ( D2HDM{xt,yt) 



47r \S2HDM{Xt,yt) 



J5 



with 



(70) 



(71) 



S2HDM{xt,yt) = [Sa{xt) + SwH{xt,yt) + SHHixt,yt)] , 
D2HDM = DsAiixt) + Dnixt^yt), 

here the high energy matching scale fj, = Mw is chosen. The functions DsM{xt) and Dh (xt , yt) of the eq. 
contain the SM and new physics parts of the NLO QCD corrections to the mass difference AMb^ [85j . 



DsAiixt) = Cf [L(i''5*^)(.Tt) + SSoixt) 



Ca 



'L^'-'''\xt) + hS,{xt) 



Dnixt^yt) = Cf L^^'"\xt,yt) + 3 {SwH{xt,yt) + SHH{xt,yt)) 



Ca 



L^^'"\xt,yt) + 5{SwHixt,yt) + SHH{xt,yt)) 



(72) 
(73) 

(74) 
(75) 
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where Cp = 4/3 and Ca = 1/3 for SU (2i)c. The function So{xt) includes the dominant top-box contribution in the 
SM and has been given in Refs. (sol l85l|. The functions SwH{xt,yt) and SHnixt^yt) incorporate the new physics 
contributions from the box diagrams with one or two charged Higgs bosons involved [85[ , 



SwH{xt-,yt) 



SHH{Vt 



^3*3^31 



VtbVtd 



ytxt 



{2xt - 8yt)ln{yt) ^ 

i^-ytyiyt" xt) {1~ xt)'^{yt~ xt) 



6xtln{xt) 



8 - 2xt 



a-yt){i 



^3*3^31 



VtbVtd 



yt xt 



Xt) 

1 



yt 



'2ytln[yt] 

{i^ytr^ {i-yty 



(76) 
(77) 
(78) 



with 



^33^31 



VtbVtd 



Y 



m 

fm, 
mt 



Vcb\f{Y) 
VtbJ V2 



X23 



Y 



V 



Vtb) 72 



X23 



A 



M 

Vcd\f{Y) 
Vtd) ^/2 



X33 



X23 



Vtd) ^/2 



(79) 



where we have used ec 
are not considered 124 



in order to compare with previous results in the literature, where the non-diagonal terms 
so that one can recover those results when X23 = 

0. Finally, the function L**'^) {i = 1,8) 



describes the charged Higgs contribution [8f 



L'''^"\xt,yt) 



^3*3^31 



VtbVtd 



WH^'\xt,yt) +2 



^3*3^31 



VtbVtd 



<i>H^'Hxt,yt) + 



^3*3^31 



VtbVtd 



HH('Hyt) 



(80) 



The explicit expressions of the rather complicated functions WH'^'^^xt,yt), ^H'^^^xt,yt) and HH'--^\yt) can be found 
iuRef. m. 

Following the analysis of Ref. [s^l and considering the areas allowed by the measured AMb^ value within a 2a 
error, when we have a light charged Higgs (80 GeV < toh± < 200 GeV) one can extract the limit 



^3*3^31 



VtbVtd 



< 0.25, 



(81) 



which is consistent with the bounds obtained from B — > Xs"f in the previous subsection. We present then Fig. |S1 
which is the allowed region in the plane [X33, X23] the cases F << 1 (left panel), Y = 1 (center panel) and Y = 10 
(right panel). One can get the bound —0.75 < X23 — —0.15 for X33 = 1, and 0.4 < X23 ^ 0-9 for X33 = —1, both when 
y << 1. This result reduces a lot the intervals of Tab. |lT]for the cases presented. 



G. 



bb 



Other stringent bounds on |x33| come from radiative corrections to t he p rocess Z — > bb, especially to the Z decay 
fraction into bb (Rb). Following the formulas presented in Refs. [6ll. |62|. [7a [86j . Rb is parametrized as: 



Rb — 



r{z bb) 

r{Z hadrons) 



= 1 + ^ 



where 



h = 

Kb = 



{9b ~ 9^? + {9b 



QCD 



q^b,t 



Kb 

h 



(82) 



(83) 



with C^^^ = 1.0086, which is a factor that includes QCD corrections, and 



-L,R 



-L,B. 
3zbb 



Sg 



LM 



(84) 



where g^'f^j^ are the tree level couplings and Sg^'^ are the radiative corrections that include the contributions of new 



physics. In models with two doublets Sg^ and Sg^ have contributions from loops involving all the Higgses {H^ , H^, 
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h° and Then, following the calculation of Ref. [t^, we obtain that the dominant contributions for Sg^'^ come 
from the charged Higgs boson and are given by 



167r2 
R 



mt 

RlogR 



Y 



R-1 (i? - 1) 



V2 

2 • MR) } , 



167r2 
R 



RlogR 



" ^^^^ 



+ /2(i?) , 



(85) 



(86) 



R-l (i?-l)2 

where R = ml/m^^^ and the function ,f2{R) governing the NLO corrections is given in Refs. |72| . Again, when 

X33 



— 0, one gets the case for the 2HDM-II. According to the measured value of i?t [87|, 



Rb = 0.21629 ±0.00066 



(87) 



we can get the experimental constraints for 5Rb = |0.00066|. Then, from eqs. (|85 |) - ([86l) we obtain bounds for Ixggj and 
X, Y . From Ref. [HH we can also use the combined limit from leptonic r decays and the global fit to (scmi-)leptonic 
decays, which is given by: 



\Y33Z3 



< 0.005, |X33|<50 



(88) 



However, this constraint is already contained in the 6 — >■ 57 ones, which are more restrictive, so that this last result 
does not modify those obtained in the previous section. 



IV. LIGHT CHARGED HIGGS BOSON PHENOMENOLOGY 



Now we discuss some simple yet interesting phenomenology emerging in the 2HDM-III after all aforementioned 
constraints are taken into account, in particular the decays of a light charged Higgs boson (i.e., with mass below 
the top quark one). The expressions for the charged Higgs boson partial decay width i7+ uidj to massive quark- 
antiquark pairs are of the form 



3.9^ 



327rAf,i 
1 
2 



-Ai/2(r 



(5„ 



Pi,) 



(89) 



where A is the usual kinematic factor A(a, b,c) = {a — b — c)^ — 4bc. When we replace Xwd ^ 0, the formulas of the 
decay widths become those of the 2HDM with NFC: see, e.g., Refs. [1, Q. For a scenario with light charged Higgs 
bosons, for which we can then neglect decays into heavy fermions, the expression for the partial widths of the hadronic 
decay modes of if='='s is reduced to: 



3GFmH±{ml\X,,\^ + ml\Y,,\^) 



(90) 



For the case of leptons, one has instead 



l^^^) = 



GFmH±mf,\ 



(91) 



Again, when the parameters xjj ^"^^ absent, we recover the results of the MHDM/A2HDM. In r(iJ± u^dj) the 
running quark masses should be evaluated at the scale of mjj± and there are QCD vertex corrections which multiply 
the above partial widths by (1 + 17as/(37r)). In the 2HDM the parameter tan/3 determines the magnitude of the 
partial widths. The BRs are well known and for the case of interest, mff± < mt, one finds that the dominant decay 
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channel is cither — >■ cs or tv, depending on the value of ta.n/3. In the 2HDM-I the BRs are independent 

of tan/3 instead and the BK{H^ tv) is about twice the BR(_ff^ cs). 

The magnitude of BR(iJ='= — > cb) is always less than a few percent in three (Type I, II and lepton-specific) of the four 
versions of the 2IIDM with NFC, since the decay rate is suppressed by the small CKM matrix element Vcb 14s)- In 
contrast, a sizeable BR(i7='= — >■ cb) can be obtained in the so called 'flipped' 2IIDM [2| for tan/? > 3. This possibility 
was not stated explicitly in (H though, where the flipped 2IIDM was introduced. The first explicit mention of a large 
BR(iJ='= — >■ cb) seems to have been in Q and a quantitative study followed soon afterwards in As discussed in 
j4l| though, the condition mfj± < rrit in the flipped 2IIDM would require additional new physics in order to avoid the 
constraints on m^± from b sj, while this is not the case in the MHDM/A2HDM. Otherwise, according to the most 
restrictive flavor constraints obtained in section II, the interaction of charged Higgs bosons with the fermions in the 
2IIDM-III with a four-zero texture in the Yukawa matrices can be written in the same way as in the MHDM/A2IIDM 
plus an small deviation, namely: 

2HDM-III MHDM/A2HDM . ,„„\ 
9H±ua, ^9H±u,d, +^9- (92) 

Therefore the phenomenology of states of the 2IIDM-III is very close to than in the MIIDM/A2IIDM, although 
not necessarily the same. We will show in particular peculiar effects induced by the off-diagonal terms of the Yukawa 
texture studied here in processes considered recently in (4]| . 



A. The dominance of the BR(JJ* — > cb) 



A distinctive signal of a state from the 2HDM-III for < mt — mi, would be a sizable BR for — >■ cb. For 

TOjy± < nit — mf,, the scenario of |X| >> |Z| in a 2HDM-III gives rise to a "leptophobic" with BR(i?='^ — )■ 
cs)+BR{H^ cb) ^ 100%, as the BR{H^ tv) is negligible (<< 1%). The other decays of H"^ to quarks are 
subdominant, with the BR(7J^ — )■ us) ~ 1% always and the BR(iJ^ t*b) becoming sizeable only for m^f± ~ mt, as 
can be seen in the numerical analysis in [l^] (in the flipped 2IIDM). Incidentally, note that the case of|X| >> 
is obtained in the flipped 2IIDM for tan/3 > 3, because |X| = tan/3 = 1/\Y\ = l/\Z\ in this scenario. 

Conversely, one can see that the configuration Y >> X,Z (this imply that Yij >> Xij,Zij, see eqs. (p4 l) -([35 |) ) 
is very interesting, because the decay —J- cb is now dominant. In order to show this situation, we calculate the 
dominant terms mcl23, mcl22 of the width T{H^ — > cb, cs), respectively, which are given by: 

/ f{Y) \ fCY) 
mcYcb = TOc>23 Vcbmc\Y ^X22j - Vtb—^^/rrHr^X23 

= VcbmcX!^2 + ^thVWn^A^3, (93) 

/ f{Y) \ f{Y) 
mcYcs = TOc>22 = Vcsmc \Y yf~^22 j " Vts—^\/mtmcX23 

= VcsmcX22 + Vtss/mm'cX23- (94) 



As Y is large and ,f{Y) = ^/l + Y^ ~ Y, then the term yY — ^-^X22j could be absent or small, when Xij = 0(1). 

Besides, the last term is very large because cx yjm^m'c, given that toj = 173 GeV, so that in the end this is the 
dominant contribution. Therefore, we can compute the ratio of two dominant decays, namely, BR(i7^ — > cb) and 
BR(Jf± cs), which is given as follows: 



BR(iJ± ^ ch) \Yt 



|2 



- BR(V ^ cs) ^ 

In this case BR(i7^ — > cb) ^ 100% (for m//± < mt — mi,, of course) so that to verify this prediction would really be 
the hallmark signal of the 2HDM-III. Therefore, we can see that the non-diagonal term X23 (o^' ''^23) cannot be omitted 
and this is an important result signalling new physics even beyond the standards 2HDMs. In a similar s piri t, one can 
study other interesting channels, both in decay and production and both at tree and one-loop level [2^ l23l I28l-l30| . 
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Another case is when X » Y ,Z, here, we get that the dominants terms 

/ f{X) \ f{x) 

nibXcb mbX23 = Vcbrrib f X - ' ^ X33 j - Vcs^-^Vm/rn^xis 

VcbrribX^s + Vcsy/mTmXis, (96) 



msX22 VcsTTls I X 



^xt2)-V,s^V^sXis 



= VcsmsX22 + Vcb\/msmb\l2- (97) 
In this case there are two possibihties. If x = 0(1) and positive then — ^^^^33^ is small and 

Rsb-^^. (98) 



Here, the BR(ff ^ — )■ ch) becomes large, again, this case too could be another exotic scenario of the 2HDM-III. The 
other possibility is when x = and negative, then 

which is very similar to the cases studied recently in [4lj. In summary, one can see two possibilities to study the 
BK{H^ — > cb): firstly, the scenarios given in eqs. ([55)1 and which are peculiar to the 2IIDM-III; secondly, the 
scenario very close to the MIIDM/A2HDM expressed in eq. (|99p . 

The CKM matrix elements are well measured, with Vcb ^ 0.04 (a direct measurement), Vcs ~ 0.97 (from the 
assumption that the CKM matrix is unitary), Vts ^ 0.04 and Vtb ^ 0.999 (from direct determination without 
assuming unitarity, as is possible from the single top-quark-production cross section) (87| . With this information, it 
is enough to determine in the cases represented by eqs. (^5]) and (|M)) the dominance of the channel — )• cb for light 
charged Higgs bosons. Otherwise, the situation given in eq. ([99)) is more delicate since the running quark masses 
and rub should be evaluated at the scale Q = to//± and this constitutes the main uncertainty in the ratio Rbs- There 
is relatively little uncertainty for rUb, with nib {Q ~ 100 GeV) ~' 3 GeV. However, there is more uncertainty in the 
value of uisi although in recent years there has been much pro gres s in lattice calculations of nis and an average of 
six distinct unquenched calculations gives rrig = 93.4 ± 1.1 MeV [S^ at the scale of Q = 2 GeV. A more conservative 
average of theses calculations, = 94 ± 3 MeV, is given in [s^. In [s^ the currently preferred range at Q 2 GeV 
is given as 80 MeV < nis < 130 MeV. Using rus = 93 MeV at the scale Q = 2 GeV (i.e. roughly the central value of 
the lattice average [H,!!^) one obtains rUs {Q = 100 GeV) 55 MeV. Taking = 80 MeV and = 130 MeV at 
Q = 2 GeV one instead obtains uis ~ 48 MeV and iris ~ 78 MeV, respectively, at Q = 100 GeV. Smaller values of nis 
will give a larger BR(iJ^ — )■ cb), as can be seen in the situation given in eq. ([M]) . Note that the value rUs = 55 MeV is 
significantly smaller than the typical values rUs ~ 150 — >■ 200 MeV which were often used in Higgs phenomenology in 
the past two decades. In essence, we emphasize that the scenario |V| >> \Y\, \Z\ with mH± < rrit — mb has a unique 
feature in that the magnitude of nis is crucial for determining the relative probability of the two dominant decay 
channels of the state. Further, this is not the case for most other non-minimal Higgs sectors with states that 
are commonly studied in the literature. 

In [l0| the magnitude of BR{H^ -> cb) in the MHDM was studied in the plane of |X| and \Y\, for \Z\ = and 
0.5, taking nis = 0.18 GeV and nib = 5 GeV. With these quark masses the maximum value for the above ratio is 
Rbs = 1.23, which corresponds to BR(i7^ — > cb) ~ 55%. However, the values of nis = 0.18 GeV and nib = 5 GeV 
are no longer realistic (as it was subsequently noted in [i^l) and two recent papers 0,11^ have in fact updated the 
magnitude of Rbs (in the fiipped 2HDM) using realistic running quark masses at the scale Q = niu±. In [l^, it 
appears that nis = 0.080 GeV at the scale Q — m^i was used, which gives BR(iJ^ — >■ cb) ~ 70%, in agreement 
with the results presented in [4l| . In nis = 0.077 GeV at the scale Q = niu± was used, with a maximum value 
for BR(iJ^ cb) of ~ 70%. We note that none of these papers used the precise average wig = 93.4 ± 1.1 MeV [s^ 
of the lattice calculations, which gives nis ~ 55 MeV at the scale of Q = mff± . This smaller value of nis leads to a 
maximum value of BR(i7* — > cb) which is larger than that given in 0, [13, 1121: as discussed below. 

We now study the magnitude of cb as a. function of the couplings X, Y, Z. In Fig. [TU] we show the numerical 

study of BR{H^ cb) in the plane [X,Y] in the 2HDM-III with \Z\ = 0.1, using x{j = 0.1, = 0.055 GeV and 
nib = 2.95 GeV at the scale of Q = m^fi = 120 GeV. With these values for the quark masses the maximum value 
is BK{H^ cb) ^ 71% (91%) considering m^o = 80 GeV {nij^a = 100 GeV), i.e., a significantly smaller (larger) 
value than BR(iJ± -> cb) - 81% [U, and both larger than BR(i?± cb) - 55% in [l3|. In Figs. [10] and HI] we also 
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FIG. 10. For the exotic scenario of eq. (|98[l we show the contours of the BK{H^ — >■ cs) (left) and BR(iJ* — s> cb) (right) 
in the plane [X, with Z = 0.1, m^o = 80 GeV, mu± = 120 GeV and Xij — 0.1. The constraint 6 — >■ S7 is shown as 
\fBl{X,Y)\ < 1.1 for Re[fBl{X,Y)] < 0, where fBl{X,Y) = is given in eq. For \ fBl(X,Y)\ < 0.7, it is when 

Re[fBl{X, Y)] < 0. We take m,(Q = m^±) = 0.055 GeV. 




display the bound from b ^ (for m/^i = 100 GeV), which is 



^33^3* 



< 1.1 for f/lY^'^ beinff real and negative 

Vtb Vtsp ° ° 



plus 



< 0.7 for 



-^^3 3^3*2 
VaVts 



being real and positive. The parameter space for BK{H^ cb) > 70% (90%) roughly 



corresponds to \X\ > 10 {\X\ > 5) and |r| < 0.35 for 



^331-3; 



VttVt, 



< 0.7. In Fig.Hniandiniwe also show BR(iJ± cs) 



The latter could be 40% (25%), however, it cannot be maximized for |y| >> \X\, \Z\, because eq. (^5]) does not allow 
it, whereas BR(i7^ — > rz^) is maximized for \Z\ >> \X\, \Y\. In Fig. [T^] we show contours of BK{H^ — > cb) and 
BR(i7± Tv) in the plane [X,Z] for m^i = 120 GeV, m^o = 100 GeV and |y| = 0.05. For this value of |r| the 
constraint from b 57 is always satisfied for the displayed range, \X\ < 20. One can see that the largest values of 
BR(i7± ^ cb) arise for |Z| < 2. 

Prospects for t — )- H^b with — > cb at the LHC have been reviewed lately in Ref. [iJl for the case of the 
MHDM/A2IIDM, from which many results can however readily be adapted to our current studies. Things go as 
follow. The case of mH± < rrit — rrih with a large BR(7J^ — >■ cs) can be tested in the decays of the top quark via 
t H^b and was studied first in 0, Is^l- Innumerable studies in this direction followed suit, far too many in fact 
for being listed here. Also, we are concerned here primarily with — >■ cb decays. The first discussion of t — >■ H^b 
followed by — >■ cb was given in [s^. Recently, t H^b with — > cb has been studied in the context of the 
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flipped 2HDM Sm and the 2HDM without NFC H (other than in (U). 

Light charged Higgs bosons {H^} are being searched for (or simulated) in the decays of top quarks (t H^b) at 
the Tevatron and at the LHC in a variety of modes [s^ - IsgI . [9lll92j. In particular, separate searches are being carried 
out for the decay channels — ^ cs and — > tv, with comparable sensitivity to the mass and fermionic couplings 
of . The searches for — > cs in [ssf and [36| look for a peak at m//± in the dijet invariant mass distribution, 
with the assumption that neither of the quarks is a 6-quark. 

In the MHDM/A2HDM reahsations the BR(i7± cb) can be as large as 80% for an that is light enough to be 
generated in the t — >■ H^b decay. This should be contrasted to the case of the state belonging to more standard 
2HDMs for which a large BR{H^ — > cb) is certainly possible but one expects mjy± > rrit in order to comply with the 
measured value of 6 — > sj, which is not the case for the 2HDM-III. The latter scenario is therefore on the same footing 
as the MIIDM/A2HDM cases studied in [4l|- Herein, in short, it was suggested that a dedicated search for t H^b 
and — > cb would probe values of the fermionic couplings of which are currently not excluded if one required a 
6-tag in one of the jets originating from , thus affording in turn sensitivity to smaller values of BR(t — > H^b) than 
those obtained to date (which use un-flavored jet samples). Therefore, a dedicated search for t H^b and cb 
at the Tevatron and LHC would be a well-motivated and at the same time simple extension of ongoing searches for 
t -> H^b with decay -> cs. 



B. The decay -> AW* for m^o < m^± 

The above discussion has assumed that cannot decay into other (pseudo)scalars. We now briefly discuss the 
impact of the decay channel — >■ ^"^VF*, which has been studied in the 2HDM-II in [o^l and in other 2HDMs 
with small |X|, |y| and \Z\ in [iOl, in the light of direct searches at LEP (assuming A° — >• bb) performed in (ol ]. 
In a general non-SUSY 2HDM the masses of the scalars can be taken as free parameters. This is in contrast to 
the MSSM in which one expects ~ m^o in most of the parameter space. The scenarios of m^o < mjj± and 

m^o > mjj± are both possible in a 2HDM, but large mass splittings among the scalars lead to sizeable contributions 
to EW precision observables [9^, which are parametrized by, e.g., the S, T and U parameters [9^. The case of exact 
degeneracy (m^o = 7ti^o = rnH±) leads to values of 5, T and U which are almost identical to those of the SM. A 
recent analysis in a generic 2HDM [97| sets m^/o = m^o, sin(/3 — a) = 1 and studies the maximum value of the mass 
splitting Am = m^o — m//± (for earlier studies see (sMl). For m^o = 100 GeV the range — 70GeV < Am < 20GeV 
is allowed, which corresponds to 80 GeV < m//± < 170 GeV. For to^o = 150 GeV the allowed range is instead 
—70 GeV < Am < 70 GeV, which corresponds to 80 GeV < < 220 GeV. Consequently, sizeable mass splittings 

(of either sign) of the scalars are possible. Analogous studies in a MHDM have been performed in [o^, with similar 
conclusions. 

If m^o < mu± then the decay channel — >■ ^''Vl/^* can compete with the above decays of to fermions, 
because the coupling H^A'^W is not suppressed by any small parameter. In Fig. [13] (left) we show contours of 
BR(H± ^ A°W*) in the plane [X,Y] with \Z\ = 0.1, m.^o = 125 GeV and tojj± = 150 GeV whereas in Fig. [H] 
(right) we present the case m^o = 80 GeV and m^± — 120 GeV. The contours are presented in the parameter 
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MA= SO GcV 
IZI=0.1 



10% 
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FIG. 13. We show the contours of the BR{H^ AW) for m^^i = 150 GeV, m^o = 125 GeV (left) and m^± = 120 GeV, 
m^o = 80 GeV (right) in the plane [X, Y] with Z = 0.1, X >> Y and x{j = 0-1 (^'U positive). The constraint 6 — > 57 is shown 
as \ fBl{X,Y)\ < 1.1 for Re[fBl{X,Y)] < (red-dashed), and \ fBl{X,Y)\ < 0.7 is when Re[fBl{X,Y)] < (blue-dashed). 
We take m^Q = m^±) = 0.055 GeV. 



space of interest (i.e., |y| < 1 and \X\ >> 1) because the contribution of the term m^|y22p to the decay widths of 
to fermions could decrease. Comparing, e.g.. Fig. [13] (right) to Fig. [10] (right) one can see that for \X\ ^ i a 
BR(i/± -> A°W*) ~ 30% is competitive with a BR(i7± cb) - 65%. For smaller m^o (e.g., less than 80 GeV) the 
contour of the BR(ff± ^ A°W*)=50% would move to higher values of \X\. Since the dominant decay of the A'^ is 
expected to be A^ — > bb, the detection prospects in this channel should also be promising because there would be 
more b quarks from t — > H^b, — )- A°W* , A^ — >• bb than from t — > H'^b with — )■ c6. We note that there has 
been a search by the Tevatron for t he channel t — >■ H^b, — >■ A°VF*, A° — >■ r+r" [100| . for the case of to^o < 2mb 
where AP — > bb is not possible |lOl| . 



At present there is much speculation about an excess of events around a mass of 125 GeV in the search for the 
SM Higgs boson [4^. An interpretation of these events as originating from the process gg ^ A^ ^ 77 has been 
suggested in |l02j . In Fig. [T3] (right) we set m^o = 125 GeV and mfj± ~ 150 GeV for this reason. Since the 
mass splitting between and A^ is less than in Fig. [13] (left), the contours move to lower values of but a 
BR{H^ A^W*) = 4% is still possible for \X\ < 2. We note that if the excess of events at 125 GeV is attributed 
to a SM-like Higgs, then in the context of a 2HDM, a candidate would be the lightest CP-even Higgs /i° with a 
coupling to vector bosons of SM strength (recent studies of this possibility can be found in |l03{ ). This scenario would 
correspond to sin(/3 — a) ^ 1 in a 2HDM and therefore the couphng H^h°W (with a magnitude cos(/3 — a) in a 
2HDM) would be close to zero. Hence the decay h°W* would be suppressed by this small coupling, as well as 

by the high virtuality of W* . Several recent studies |104 Il05l | fit the current data in all the Higgs search channels 
to the case of a neutral Higgs boson with arbitrary couplings. A SM-like Higgs boson gives a good fit to the data, 
although a slight preference for non-SM like couplings is emphasized in [105[. If the excess of events at 125 GeV 
turns out to be genuine and is well described by a non-SM like Higgs boson of a 2HDM-III with a value of sin(/3 — a) 
which is significantly less than unity, then the BR(iJ^ — > h^W*) could be sizeable, with a magnitude given by Fig.[T5l 
(right) after scaling by cos^(/3 — a). 



NUMERICAL RESULTS 



A. Indirect production 



We now quantify the magnitude of — >■ cb events produced in the decays of t quarks, and compare this with the 
expected sensitivity at the LHC. For the partial decay widths of t — W^b and t H^b we use the LO expressions 
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(with \Vtb\ = 1) as follows: 



r{t ^ WH) = + 2A/^][l - A/^/™?]2 (100) 



Tit ^ HH) = |£^[„.?|y33P + ml\XMl m\^lm\\\ (101) 

The multiplicat ive ( vertex) QCD corrections to both t — > W^h and t — > R^h essentially cancel out in the ratio 
of partial widths jl06l |. In the phase-space function of both decays we neglect m;,, and in the terms mjlYsap and 
TO^|X33p we use nit = 175 GeV and TOf, evaluated at the scale m^± (i.e., mf, ^ 2.95 GeV). 

In Fig. [14] we show contours of the sum of 

BR(t ^ HH) X [BR(i?± ^ cs) + BR(F± ^ c6)] (102) 

in the plane of F] for m//± = 120 GeV and m//± = 80 GeV, setting \Z\ = 0.1. The cross section in eq. (|102p 
is the signature to which the current search strategy at the Tevatron and the LHC is sensitive, i.e., one 6-tag (LHC 
[3^) or two 6-tags (Tevatron [s^) are applied to the jets originating from the tt decay, but no 6-tag is applied to the 
jets originating from H^. For the case of \SK{H^ — >■ cs)-|-BR(i7^ — > c6)]=100% the current experimental limits for 
m^± = 120 GeV are BR(t ^ HH) < 0.14 from ATLAS with 0.035 fb'^ [s^, BR{t -> HH) < 0.12 from CDF with 
2.2 fb-i Hi, and BR(< -> iJ±6) < 0.22 from DO with 1 fb'^ [H. In Fig. HI] (left) for m^i =120 GeV these upper 
limits would exclude the parameter space of |X| > 50 and small |y| which is not excluded by the constraint from 
b 57. For the mass region 80 GeV < m^i < 90 GeV there is only a limit from the DO search in [sj], which gives 
BK{t — )■ H^b) < 0.21. From Fig. [14] (right), for m^± = 80 GeV, one can see that this limit excludes the parameter 
space with \X\ > 45 and small 

In Fig. [T4]we show contours of 1%, which might be reachable in the 8 TeV run of the LHC. Simulations by ATLAS 
(with Vs = 7 TeV) for ^ cs ^ have shown that the LHC should be able to probe values of BR{t HH) > 0.05 
with 1 fb~^ for m}j± > 110 GeV, with the greatest sensitivity being around to^± = 130 GeV. For the operation 
with y/s = 8 TeV and an anticipated integrated luminosity of 15 fb~^ one expects increased sensitivity (e.g. BR(< — >■ 
H^b) > 0.01 for mH± > 110 GeV), although the region 80 GeV < mH± < 90 GeV might remain difficult to probe 
with the strategy of reconstructing the jets from H^. An alternative way to probe the region 80 GeV < rn^i < 90 
GeV is to use the search strategy by DO in [3^ . and presumably the LHC could improve on the Tevatron limit on 
BR(t — ^ H^b) < 0.21 for this narrow mass region. From Fig. [14] (left) (for mH± = 120 GeV) one can sec that 
the region of |F| > 0.32 and \X\ < 14, which is not excluded by 6 -> sj, would be probed if sensitivities to the 
BR(t — >■ H^b) > 0.01 were achieved. However, a large part of the region roughly corresponding to \Y\ < 0.32 and 
\X\ < 20 (which is also not excluded by & — >■ 57) would require a sensitivity to BR(t — > H^b) < 0.01 in order to be 
probed with the current search strategy for t — )- H^b and this is probably unlikely in the 8 TeV run of the LHC. 

Increased sensitivity to the plane of [X, Y] can be achieved by requiring a &-tag on the jets which originate from 
the decay of H^. In Figs. ITS] and ITCl for m^± = 80 GeV and m/^± — 120 GeV, respectively, we show contours of 

BR{t hH) X BB.{H^ -> cb) . (103) 

With the extra 6-tag, as advocated previously (see eq. (12) in [4l|), the sensitivity should reach BR(t H^b) x 
BR{H^ cb) > 0.5%, and perhaps as low as 0.2%. In the latter case, one can see from Figs. [TS] (left) and [TBI fright) 
that a large part of the regions of \X\ < 6 (for m//± = 120 GeV) and \X\ < 4 (for TOjy± = 80 GcV) could be probed, 
even for \Y\ < 0.2. Therefore, there would be sensitivity to a sizeable region of the parameter space of [AT,!"] which 
is not excluded by 6 —> 57, a result which is in contrast to the above case where no 6-tag is applied to the 6-jets 
originating from . We encourage here a dedicated search for t — > H^b and — cb by the Tevatron and LHC 
collaborations. Such a search would be a well-motivated extension and application of the searches which have already 
been carried out in [lH and [s^ and would offer the possibility of increased sensitivity to the fermionic couplings and 
mass of in not only the MDHM/ A2HDM but also the 2HDM-III. 



B. Direct production 

Another possibility to produce states in our scenario is via cb (and 6c) fusion in hadron-hadron collisions. 
Since neither of the antiquarks in the initial state is a valence state, there is no intrinsic advantage in exploiting 
proton-antiproton collisions at the Tevatron, further considering the much reduced energy and luminosity available 
at the FNAL accelerator with respect to the LHC. Hence, we focus our attention onto the latter. In fact, we can 
anticipate that, amongst all the energy and luminosity stages occurred or foreseen at the CERN machine, only the 



26 





X X 
FIG. 15. Contours of BR{t H^h) x BR(H± ^ cb) in the plane [X, Y] with \Z\ = 0.1 for m^i = 80 GeV. The constraint 
6-^ S7 is shown as \ fBl{X,Y)\ < 1.1 for Re[fBl{X,Y)] < (red-dashed), and |/_B1(X,F)| < 0.7 is when Re[fBl{X,Y)] < 
(blue-dashed). We take ms{Q = mn±) = 0.055 GeV and show the range < |y| < 0.8 (left panel) and < |y| < 0.3 (right 
panel) . 



combinations y/s — 8, 14 TeV and standard instantaneous luminosity of order 10"^^ cm^^ s^^ or higher are of relevance 
here. Furthermore, in the QCD pohuted environment of the LHC, it is clear that one ought to attempt extracting a 
leptonic decay of our hght charged Higgs boson. In the light of the results presented in the previous (sub)sections, 
the best option is afforded by decays into riy pairs, eventually yielding an electron/muon (generically denoted 
by I) and missing (transverse) energy. In this case, the background is essentially due to the charged Drell-Yan (DY) 
channel giving riy pairs (i.e., via production and decay), which is in fact irreducible. 

Figs. [TTHTSI displav the differential distribution for the signal S (i.e. cb t^v + c.c.) and the background 

B (i.e., c6 — )■ W^T^v -\- c.c), the former for three mass choices {mfj± — 80, 120 and 160 GeVjl. As the invariant 
mass of the final state is not reconstructible, given the missing longitudinal momcmtum, wc plot the transverse 

mass Mt = \J {Ef + E'^^^^^)'^ — {pf + pf^iss)"^ ~ (pf +-Pmiss)^' where E'^ represents missing energy /momentum (as we 
consider the electron and muon massless) in the transverse plane and Px,y are the two components therein (assuming 
that the proton beams are directed along the z axis). Clearly, the backgound dominance over the signal is evident 
whenever m[j± sa M^y± . However, the larger the charged Higgs boson mass with respect to the gauge boson one, the 
more important the signal becomes relatively to the background. 

In order to establish whether it is possible to extract the direct signal, or indeed constrain it, we show in 
Figs. fT9l-Figs.[20l the significance of the signal, defined as S/VB, as a function of the collider integrated luminosity L, 



Notice that interference effects are negligible between the two, owing to the very narrow width of the charged Higgs boson. 
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FIG. 16. Contours of BR{t H^b) x BR(_ff± cb) in the plane [X, Y] with \Z\ = 0.1 for m^i = 120 GeV. The constraint 
b ^ S7 is shown as \ fBl{X, Y)\ < 1.1 for Re[fBl{X, Y)] < (red-dashed), and \ fBl{X, Y)\ < 0.7 is when Re[fBl{X, Y)] < 
(blue-dashed). We take ms(Q = m^±) = 0.055 GeV and show the range < |y| < 0.8 (left panel) and < |y| < 0.3 (right 
panel) . 
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FIG. 17. Differential distributions in transverse mass Mt for signal and background (the former for three mass values) in 
logaritmic (top) and linear (bottom) scale. Here, ^/s = 8 TeV. 



where S = L a{cb t^v + c.c.) and B = L a{cb — > W^t^v + c.c). Here, we have restricted both the signal 

and background yield to the mass regions \mH± — Mt\ < 10 GeV, where, again, we have taken mH± = 80, 120 and 
160 GeV. The t decay rates into electron/muons are included and we assumed 90% efficiency in e//x-identification. It 
is clear that exclusion (significance equal to 2), evidence (significance equal to 3) and discovery (significance equal to 
5) can all be attained at both the 8 and 14 TeV energy stages, for accessible luminosity samples, so long that the 
mass is significantly larger than the one. In both machine configurations, corresponding event rates are always 
substantial. 



VI. CONCLUSIONS 



Light charged Higgs bosons (i?^) are being searched for in the decays of top quarks {t — >■ H^h) at the Tevatron 
and LHC. Separate searches are being carried out for the decay channels — >■ cs and — >■ tv^ with comparable 
sensitivity to the mass and fermionic couplings of H^. The searches for — )- cs in [sH] and [s^ look for a peak at 
mfj± in the dijet invariant mass distribution, with the assumption that neither of the quarks is a 6 quark. 

In some models with two or more Higgs doublets (the 2HDM and a MHDM with three or more scalar doublets) the 
BR for iJ± cb can be as large as 80%. Here, in our model (2HDM-III), BR(7?± cb) could be as large as 90%. 
Moreover, such a could be light enough to be produced via t H^b as well as respect the stringent constraints 
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FIG. 18. The same as Fig. [T71 but for ^=U TeV. 
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FIG. 19. Number of events (top) and significance (bottom) of the signal (for three mass values) as a function of the 
luminosity. Here, -^5 = 8 TeV. The horizontal solid(dashed) [dotted] line corresponds to a significance of 5(3) [2]. 



from 6 — >■ S7 on both mu± and the fermionic couphngs of . This is in contrast to the state in other 2HDMs 
for which a large BR for H'^ — >■ c6 is possible but for which one expects ra^i^ > rrit in order to comply with the 
measured value of 6 — > 57. Along the lines indicated by previous literature, in the context of the 2HDM-III with a 
four-zero Yukawa texture, we suggested that a dedicated search for t — ?> H^h and cb would probe values of 

the fermionic couplings of that are currently not testable. Such a search would require a 6-tag of one of the jets 
originating from and would further allow one to reach a higher sensitivity to a smaller value of the BR(t — > H^b) 
than that obtained in the ongoing searches, which currently do not make use of this additional 6-tag. 

Finally, we have shown that a state of the 2HDM-III can also be produced directy from cb fusion, followed 
by a decay into tv pairs, assuming semi-leptonic decays (to contrast the otherwise overwhelming QCD background), 
albeit only at the LHC, which should have sensitivity to it already at 8 TeV and, indeed, full coverage at 14 TeV, so 
long that TTiHi is 120 GeV or above. 

We have reached these conclusions over the 2HDM-III parameter space that we have shown to survive a long list 
of constraints emerging from B-physics, namely: fi — e universality in r decays, several Icptonic _B-decays {B — >■ ri^, 
D ^ fiv and Dg Iv), the semi-leptonic transition B — >■ Dtv, plus B — > X^j, including — B^ mixing and the 
radiative decay Z — > bb. 

The outlook is therefore clear. Depending on the search channel, both the Tevatron and the LHC have the potential 
to constrain or else discover the 2HDHM-III supplemented by a four-zero Yukawa texture including non-vanishing 
off-diagonal terms in the Yukawa matrices. We are now calling on our experimental colleagues to achieve this, as the 
search strategies recommended here are easily implementable and pursuablc. 
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FIG. 20. The same as Fig. \W\ but for y/s = U TeV. 
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